#!/usr/bin/python

from ceorModule import *
import time

######################################
################ MAIN ################
## Is a 2-D box with lenght xL yL
if __name__ == "__main__":
  # And God says...
  eps = 1.0
  sigma = 1.0
  delta_r_ini = 2.5*sigma
  dt = 0.001
  gas = Gas2D(iniFile="N_000002_Lx_100_Ly_10_t_0000.dat")
  tic = time.time()
  gas.initPos()
  gas.arrayPos()
  gas.arrayVel()
#  for pi in gas.particle:
#    pi.vy_new = 0.0
  toc = time.time()
  print toc-tic
#  gas.plot()
  print gas.particle[0].x_new, gas.particle[1].x_new
  
  gas.particle[0].x_new = 1.0
  gas.particle[1].x_new = 2.5
  gas.calculateForces(eps, sigma)
  
  print gas.particle[0].fx_new
  print gas.particle[1].fx_new
#  gas.loadtxt()
